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EVALUATION OF THE THIN DEFORMABLE ACTIVE OPTICS MIRROR CONCEPT 

By Hugh J. Robertson 
The Perkin-Elmer Corporation 


1.0 INTRODUCTION 

1.1 OBJECTIVES 

The Active Optics concept using a thin deformable mirror has been 

successfully demonstrated using a 30" diameter, 1/2" thick mirror (figure 1-1) 

* 

and a 61 point matrix of forces for alignment . Many of the problems associ- 
ated with the design, fabrication, and launch of large aperture diffraction- 
limited astronomical telescopes have been resolved and experimental data cre- 
ated that can provide accurate predictions of performance in orbit. 

The work that has been performed under this contract was designed 
to further evaluate the present system to determine its parameters more com- 
pletely in order to facilitate the optimum design for larger active telescopes. 
Specifically the evaluation included the following tasks as design objectives: 

1. Determination of the maximum degree of control possible as a 
function of actuator spacing. (The range of radius of curvatures that the 
mirror can assume with a given number of actuators.) 

2. Experimental determinations of mirror influence coefficients. 

3. Determination of the effects of thermal gradients on the thin 
mirror and the ability of the control system to accommodate these effects. 

4. Determination, implementation, and evaluation of an actuator 
arrangement designed to take advantage of natural modes of deflection of the 
mirror. In particular it was desired to determine the optimum positioning of 
actuators and distribution of actuator forces in order to minimize the number 
of control channels required and to decouple their action so as to optimize 
system response. 

* “ 
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Figure 1- 


30-Inch Diameter, Thin, Deformable Mirror 
and Actuator Assembly 
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5. Determination of the response of the control 
system to step error disturbances with and without 
feedforwards for each actuator arrangement. 

6. Determination of the degree to which localization 
is achieved using the force patterns derived analy- 
tically under NASA Contract NAS 1-7103. 

1.2 SUMMARY OF RESULTS 

1.2.1 Range of Mirror Control 

The range of radius of curvatures that the active thin mirror 
could assume was investigated and the figure error was plotted against change 
in radius of curvature using 4, 16, and 58 actuators. The mirror figure was 
observed to deteriorate exponentially as the radius varied from the minimum 
error alignment position, at a rate inversely proportional to the number of 
actuators used. The 58 actuators with a maximum force capability of ±2 pounds 
were able to change the radius of curvature over a range of approximately 1/2 
inch. The effects of the changes of curvature on the mirror figure and its 
ability to image a pinhole source were documented using interferograms , pinhole 
photographs and pinhole image profile scans. 

1.2.2 Experi menta l De ter mination of Influence Coefficients 

To obtain accurate experimental data on the 30-inch thin mirror 
in order to provide more precise calibration of analytical determination tech- 
niques, the mirror influence coefficients were measured. The mirror influence 
coefficients were determined experimentally by measuring the displacement of 
the stressed mirror at various points, from interferograms. The necessary 
data was acquired and reduced to give a full set of the individual rows of the 
influence coefficient matrix. 

1.2.3 Thermal Measurements 

The thermal effects induced by the space environment and the 
resulting optical performance of a large aperture diffraction limited system 
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is a major problem area that bears more detailed analytical and experimental 
investigation. Once the anticipated temperature and temperature changes have 
been established, the resulting effect on performance can be predicted and 
proper open or closed loop compensation can be designed into the system and/or 
spacecraft. As yet, however, no complete analytical solution to the thermo- 
elastic problem in large primary mirrors has been developed. Among other things, 
the solution requires the development of mathematical models of complex mirror 
structures. This represents a tedious, time-consuming and extremely complex 
analytical problem and can best be handled via suitable approximations once 
parameters can be experimentally determined. 

The thermal program was undertaken to experimentally determine 
figure degradation and focus shifts resulting from various thermal inputs to 
a thin, 30-inch diameter, 1/2-inch thick deformable mirror. The induced thermo- 
elastic deformations of a 30-inch diameter, 1/2-inch thick, fused silica pri- 
mary mirror were recorded with a phase measurement interferometer, and the ability 
of the active control system to maintain alignment in the presence of thermal 
gradients was assessed. Various thermal patterns were developed in the test 
mirror by 26 electric heater pads located at the mirror back and edge surfaces. 
The temperatures at the mirror surfaces were measured by approximately 105 
five-mil thermocouples that were attached to the mirror. The tests were con- 
ducted in a horizontal vacuum tank with the mirror supported in the Active 
Optics System to simulate the gravity-free environment experienced in space. 

The results of the thermal measurements show that localized 
gradients generated by small area heat sources on the back of the mirror pro- 
duce localized deformations. It should be possible to predict the effect of 
distributed gradients of low spatial frequency by linear superposition of the 
deformations due to localized gradients that add up to the distributed gradient. 

1.2.4 Modal Control Using Natural Modes of Mirror Vibration 

Using the modal control approach* the final figure error of a 
deformable mirror with a given actuator configuration in an active control 


system, can be decreased relative to the final error achieved by the technique 
^Reference No . 2 . 



of nulling the figure error in the vicinity of each actuator. The corollary 
to this is that the number of actuators required to reduce an initial figure 
error a specific amount to reach a given design figure error can be signifi- 
cantly less for a modal control system than for a straight nulling system. 

The technique of modal control was investigated using the natural 
modes of vibration of the mirror as an orthogonal set. Analysis of the natural 
vibrational modes of a mirror simply supported at three points yielded solutions 
that showed good agreement with modal contours experimentally observed by time 
lapse holography. The initial figure error of the 30-inch thin mirror was 
analyzed in terms of the vibrational modes, and the number and location of 
actuators required to reduce the largest modal components was investigated. 

Because of the complexity of the problem for any large number 
of actuators, a limited number of actuator and sensor points were selected to 
demonstrate the performance of the control approach. The controller for a 
four actuator and seven sensor point arrangement was tested, and the results 
established that the modal control approach performs as predicted. 

1.2.5 Control System Response to Step Error Disturbances 

The control system response was investigated using step errors 
introduced at individual actuators and step errors introduced in the axial 
spacing between mirror and figure sensor. Errors were applied to the system 
for straight parallel control, feedforward decoupling, and modal control con- 
figurations . 

1.2.6 Feedforward Localization 

The degree to which localization was achieved using the 19 
point force patterns generated under NAS 1-7103 was investigated. It was 
found that the force patterns were not being achieved due to actuator dead 
zones. Readjustment of driving point resistance values provided a very high 
degree of localization, however, by producing seven-point force patterns. 
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2.0 RANGE OF MIRROR CONTROL 


It might be desirable for some applications for an active system to 
be able to change the radius of curvature of the mirror, possibly enough to 
produce a more useful aspheric surface. As will be shown, a spherical profile 
can be maintained over a range of radii. The exact range that can be obtained 
with a given number of actuators before the mirror figure deteriorates as a 
result of a buildup of excessive ripple between the actuators was investigated. 

The performance of the mirror as it approaches the limits of control was evaluated 
primarily from interferograms and scanning of the pinhole image. 

The pinhole image scanner, shown in figure 2-1, consists of a high 
power microscope objective which enlarges the image and refocuses it on a pin- 
hole. A motor driven micrometer moves this pinhole across the image, with the 
transmitted light entering a photomultiplier tube. A chart recorder records 
the photomultiplier output. Repeated scans provide an accurate mapping of the 
image intensity profile. 

A low power objective, beamsplitter, and eyepiece is used to view 
the image by eye and allow coarse alignment of the pinhole. A field lens in 
front of the photomultiplier assures that the light entering the tube will al- 
ways fall on the same part of the tube face, eliminating the effect of non-uniform 
response of the photocathode. 

The pinhole is removable to allow viewing the whole image and also 
to permit different sized pinholes to be used. 

The figure error of the mirror is a minimum at only one axial position. 
This position was found by computer analysis of straight-fringe interferograms 
taken with the mirror aligned at different points throughout the focus range of 
the system. This was done first with the whole system of 58 actuators operating. 
After the best-focus position was established, the mirror was aligned at that 
position and the focus was again changed with only 16 actuators active and finally 
with only four actuators. 
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Figure 2-1. Scanner, Optical Schematic Diagram 
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The rms figure error observed is plotted versus the change in radius 
of curvature of the thin mirror, produced with 58, 16, and 4 actuators, in fig- 
ure 2-2. The rate of deterioration can be seen to be approximately inversely 
proportional to the number of actuators controlling, as might be expected simply 
from geometrical consideration of the spacing between actuators. 

The asymmetry of the curves is probably due to the initial figure 
error of the mirror which has a different radius in the central zone of the 
mirror than in the outer zone so that the best fits for these zones may occur 
at different axial locations of the mirror. 

The spring constant of the central actuator was increased to extend 
the range of radii available. However, the force range of peripheral actuators 
was still not sufficient to carry the change in radius to the same level of 
deterioration on both sides of the best figure location for the case of the 58 
actuator control. 

Because of the time required to process and analyze interferograms , 
another method of quickly evaluating the mirror figure was also used. This involved 
using the image scanner to record the intensity profile of a point source imaged 
by the mirror. A number of these scans, with photographs of the image and asso- 
ciated interferograms taken under the same or similar conditions, are shown in 
figures 2-3 through 2*9 f 

As can be seen in the interferograms, and the pinhole photos and 
profiles, the errors occur in hexagonal patterns as a result of the actuator 
location pattern, and a large part of the error can be attributed to the uncon- 
trolled area at the mirror edges. 

The image quality varies greatly as the mirror figure changes and, 
for small degradations (< — \ rms) , 

~ 1 " ( T )2 ( * )2 (1) 
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Figure 2-3 Pinhole Image Profile Scans at Best Mirror 
Alignment with 58 Actuators 
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Figure 2-5. Interferogram, pinhole Image, and Pinhole Image Scans for Mirror Aligned with 58 
Actuators, 0.100 u from Position of Best Alignment (Scale in Center 1/10 Scale of 
Rest of Image Scans) 





Figure 2-4. Interferogram, Pinhole Image, and Pinhole Image Scans for Best 
Mirror Alignment with 58 Actuators (Scale in Center 1/10 Scale 
of Rest of Image Scans) 






Figure 2-5. Interferogram, Pinhole Image, and Pinhole Image Scans for Mirror Aligned with 58 
Actuators, 0.100 u from Position of Best Alignment (Scale in Center 1/10 Scale of 
Rest of Image Scans) 
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where I = maximum intensity at center of image normalized to 

m 

that of a perfect image 
X = rms deformation of mirror surfaces 

Thus, changes in the mirror error can be determined rather quickly 
by monitoring the peak intensity of the po'imt source image. This assumes, 
however, that the total intensity is also known at all times. For the present 
scans,, the intensity of the source illuminating the pinhole (He-Ne laser) was 
assumed to- be constant although it appears that it did change over long periods 
of time- For future analysis a means of measuring the to-tal intensity could be 
added to provide more quantitative inf omna-tion- 
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3.0 EXPERIMENTAL DETERMINATION OF INFLUENCE COEFFICIENTS 


Influence coefficients for the 30-inch thin mirror were previously 
obtained analytically from a structural analysis program (SAMIS)*. These in- 
fluence coefficients had never been corroborated experimentally although there 
was good reason to believe that they were reasonably accurate since a comparison 
of calculated and measured influence coefficients for a scaled model had shown 
good agreement. It was now possible, however, to make measurements of the static 
deflections of the 30-inch-diameter, 1/2-inch-thick mirror including the effects 
of the mounting structure, the mounting pads, and the other actuators. 

The influence coefficients were measured from inter ferograms of the 
mirror obtained using the phase measurement interferometer. The mirror was 
first aligned using the Active Optics closed loop control system. The control 
system was turned off and then the mirror was stressed by applying a force at 
one actuator point. The displacement at all points of the mirror with respect 
to any reference point on the mirror were measured directly from the interfero- 

o 

gram, in fringes and fractions of a fringe of the 6328A light used to obtain 
the interferogram. 

Figure 3-1 shows the interferogram obtained by displacing the center 
actuator so as to change the force it applies by three pounds with respect to 
the aligned condition. The interference pattern for the interferogram shown 
was adjusted to put the three support points in approximately the same phase to 
facilitate comparison of the mirror deflections with the calculated values, 
since the calculated values assume that the support points are perfectly rigid 
and have zero displacement. This is not exactly the case in practice, and the 
support points actually show some displacement with respect to the reference. 
Different reaction support arrangements to increase the stiffness at these 
points have been investigated. 


^Reference No. 3 
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Figure 3-1. Interf erogram of Mirror with 3-Pound Force Applied 
at Center after Alignment 
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An attempt was made to obtain a measurement of the forces required 
to obtain local mirror deflections , for comparison with the calculated stiffness 
coefficient matrix. The mirror was first aligned with the Active Optics closed 
loop control system. One channel at a time was then opened and the mirror was 
driven out of alignment at the corresponding actuator while alignment at all 
other points was maintained. The change in actuator forces should give the 
forces required to obtain a local displacement and can be determined from the 
change in displacement of the actuator nuts. A linear potentiometer attached 
to each actuator provides the actuator nut displacement measurement. A set of 
forces for a local displacement at one point provides the stiffness coefficients 
for one column of the stiffness matrix. A set of forces for a local displace- 
ment of one fringe at the center point of the mirror is shown in figure 3-2 . 

It was found, however, that it was very difficult to obtain reproducible re- 
sults for any one column of the stiffness matrix. This is believed to be due 
to dead zones in the actuators due to friction and even though the dead zones 
are equivalent to less than ±1/200 wavelength in their effect on nulling at 
any one actuator, they still allow for different sets of forces and their cor- 
responding stiffness coefficients. This effect is discussed further in Section 
6 . 0 . 


Figure 3-1 shows the inter ferogram for the deflection due to a 
centrally applied force. Figure 3-3 shows typical inter ferograms at other lo- 
cations for an applied force of one pound. Figure 3-4 shows the locations at 
which the force is being applied in figure 3-3. The actuator numbering code 
used here is shown in figure 3-5. 

The displacement of the mirror measured at all the actuator locations 
gives the influence coefficients for those locations. The displacement at any 
particular point is measured with respect to a reference point by counting the 
number of fringes and fractions of a fringe between it and the reference point. 
The three mirror support points located at 120-degree intervals around the mir- 
ror edge are taken to be in the zero displacement reference plane and the inter - 
ferometer is aligned in tilt to put these three points as closely as possible 
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Figure 3-2. Stiffness Coefficients (For Local Displacement at Center of Mirror) 




Figure 3°3. Interferograms of Mirror Deformed by One Pound Force 
at Various Actuator Locations 
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Figure 3-5. Typical Influence Coefficient Reduction from Interferogram 



at equal phase in the inter ferograms . The displacement counted in fringes is 
converted to microinches by the relation (me fringe equals 12.46 microinches. 
Figure 3-5 is a set of influence coefficients measured for the mirror deforma- 
tion shown in frame 6 of figure 3-3 for a force applied at point no. 17. The 

o 

deformations are given in fringes of 6328A light and in microinches and repre- 
sent an average of measurements from three interferograms . Further averaging 
was possible in some instances due to symmetry. The averaging tends to reduce 
the effect due to errors in the initial alignment and differences in the 
actuator spring constants. The one-pound force in each case is applied by dis- 
placing the actuator spring a specific distance and since all the springs do 
not have exactly the same spring constant, there was some variation in the 
actual force applied. 

The set of influence coefficients shown in figure 3-5 represents one 
row of the influence coefficient matrix and from symmetry can represent several 
other rows of the matrix by appropriate changes of coordinates. 

A full set of influence coefficients as required to produce a complete 
matrix was obtained. Inconsistencies due to differences in the spring constants 
were compensated for, after measurement of the spring constants, in the final 
matrix rows which are shown in Figures 3-6 (a) through 3-6 (n). There are only 
14 nonredundant rows, and the actuator locations (according to the numbering 
arrangement shown in figure 5-5) to which each of the rows applies is given in 
the title of each figure. The number of the actuator at which the force is 
applied gives the first subscript for the influence coefficients in the matrix 
row. The number of the actuator location at which the mirror deflection is ob- 
served gives the second subscript identifying the influence coefficient column. 
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'•46.4 # 44.4 *37.1 *22.5 (*) 0 

\ / 

\ 52.5 . 53.0 ,46.7 J 33.3 . 
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, \ 61.8 . 62.8 Jbb.l , . 

\ / 

• 72.3 • 71.4 • • • 


\ / 



Figure 3-6(a). Measured Influence Coefficients (Microinches Per Pound) Fox- 

Force Applied at Actuator #20 (Per $ Arr. Shown in Figure 5-5) 
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Figure 3-6(b) . Measured Influence Coefficients (Microinches Per Pound)For 
Force Applied at Actuator #1, 3, 43, 51, 52, or 58 
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Figure 3-6(c). Measured Influence Coefficients (Microinches Per Pound) for 
Force Applied at Actuator 2, 44 , or 50 
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Figure 3-6(d) . Measured Influence Coefficients (Microinches Per Pound) for 
Force Applied at Actuator 4, 8, 34, 38, 53, or 57 
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Figure 3-6(e). Measured Influence Coefficients (Microinches Per Pound) for 
Force Applied at Actuator 5, 7 , 36, 42, 45, or 49 
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Figure 3-6 (f). Measured Influence Coefficient (Microinches Per Pound) for 
Force Applied at Actuator 6 , 37 , or 4-1 
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Figure 3-6(g). Measured Influence Coefficients (Microinches Per Pound) for 
Force Applied at Actuator 9 , 15, 25, 33, 54, or 56 
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Figure 3-6(h). Measured Influence Coefficients (Microinches Per Pound) for 
Force Applied at Actuator 10, 14, 26, 32, 46, or 48 
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Figure 3-6 ( i) • Measured Influence Coefficients (Microinches Per Pound) for 
Force Applied at Actuator 11, 13, 27, 31, 35, or 40 
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Figure 3-6(j). Measured Influence Coefficients (Microinches Per Pound) for 
Force Applied at Actuators 12, 28, or 30 
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Figure 3-6(k) . Measured Influence Coefficients (Microinches 
Force Applied at Actuator 16, 24, or 55 
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Figure 3-6(1). Measured Influence Coefficients (Microinches Per Pound) 
for Force Applied at Actuator 17, 23, or 47 
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Figure 3-6 (m) . Measured Influence Coefficients (Micrometers Per Pound) for 
Force Applied to Actuator 18, 22, or 39 
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Figure 3-6(n). Measured Influence Coefficients (Microinches Per Pound) for 
Force Applied at Actuators 19, 21 , or 29 
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4.0 THERMAL MEASUREMENTS 


The Active Optics system in space will be subject to thermal 
disturbances due to many heat sources, such as the sun’s irradiation on the 
primary mirror or on the side of the telescope, and electrical and mechanical 
components being heated during operation. The thermal disturbances will cause 
distortion of the mirror figure, and the function of the Active Optics system 
is to correct for distortions and maintain the performance of the diffraction- 
limited optics. 

The object of the thermal tests was to investigate the performance 
of the Active Optics system in the presence of thermal disturbances and to 
determine the response of the thin mirror to various applied thermal gradients. 
In the experiment, a given set of thermal gradient profiles was imposed on the 
back surface of the mirror, and for each of the profiles the effects of these 
gradients on the curvature of the mirror were observed and measured quantita- 
tively with the phase measurement interferometer both with and without active 
control, and the ability of the closed loop control system to maintain the 
mirror figure was assessed. In this approach the laboratory environment does 
not affect the experimental results because the experimental data for a given 
temperature distribution on the mirror are independent of the pressure of the 
test tank, the ambient temperature of the laboratory, the size and geometry 
of the test setup, and the various thermal transfer mechanisms by which the 
given temperature distribution on the mirror is finally achieved. 

The set of thermal gradient profiles that were to be imposed on 
the mirror consisted of simple and basic profiles such as a constant or linear 
or radial gradient , as shown in figure 4-1 (a), and localized thermal gradients 
produced by heat applied at small areas on the back of the mirror as shown in 
figure 4-l(b). It was assumed that most arbitrary thermal distributions can 
be considered as linear combinations of these basic profiles and that the re- 
sulting mirror deformation is the sum of the deformations due to the basic 
thermal profiles. In practice, the distributed thermal profiles shown in 
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figure 4-l(a) were difficult to obtain as uniformly as desired so that some 
smoothing was required by averaging measurements around the axes of symmetry. 

The design of the heating blanket used to produce the desired 
thermal gradients is shown in figure 4-2 and the design of the guard heater 
is shown in figure 4-3. The guard heater is located around the edge of the 
mirror to maintain distributed gradients that might otherwise be difficult to 
obtain because of heat loss to the test tank which is very close and is essen- 
tially a room temperature heat sink. The heating blanket and guard heater 
consist of many independently heated regions in which the temperature is con- 
trolled to achieve the desired temperature distribution on the mirror surface. 
The heating pads are several layers of silicon rubber that enclose heating 
elements in separate regions, and the nominal power dissipation is 1 watt per 
square inch. The temperature measurements were made by 5 mil cu/constantan 
thermocouples and thermistors, achieving an accuracy of 0.1 °C differential 
temperature . 


The thermocouples were attached to the front and back surfaces of 
the mirror by a patch of pressure sensitive aluminum tape which equalizes the 
temperature of the surrounding area of the thermocouple to the mirror surface 
temperature. The thermocouple readings were made by a Leeds-Northrop K-5 
potentiometer with a melting ice thermal bath as the reference. A test panel 
was constructed to facilitate the sampling of the thermocouples. 

Some items of interest that were considered in connection with 
the thermal test include the effects of the pressure, the radiation, and the 
profiles of the thermal gradient. 

In the actual space environment the maximum pressure will be 10 ^ 
torr for 150-Km altitude and will be lower at higher altitude (see figure 4-4). 
The effect of the presence of higher pressure near the surface of the mirror 
in the experiment is to introduce larger convection heat transfer and, thereby, 
to increase the total heat flow rate. The simplified analysis of the free 
convection heat transfer for the vertical plane indicates that 


42 






oh# 


Geometric Altitude (km) 



_ ^ 

Pressure (mb) or Density (gm ) 


Figure 4-4. Profiles to 700km of Temperature, Pressure, and Molecular 
Weight, According to U.S. Standard Atmosphere, 1962 
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where 


N NU & [<»PR »GR> 1 

n gr ■ # [p 2 ] 

N T = Nuselt's number 
NU 

N - Grashof number 
GR 

N = Prandtl number 
PR 

P = Pressure 

xf (x) and 9 (x) = sine functions of x 

Therefore , the free convection heat transfer coefficient is approximately 
proportional to the square root of the pressure (figure 4-5). As the pressure 
is decreased in the experiment, the localization of the temperature gradients 
becomes easier to maintain until the effect of the pressure becomes negligible 
compared with the radiation heat transfer through the vacuum tank wall. 

The advantage of testing the Active Optics system for a given 
temperature distribution was that the measurement data were not dependent on 
the pressure or the radiation heat transfer through the black walls of the 
test tank. 

4 . 1 THERMAL INSTRUMENTATION 

Five^mil copper/constantan thermocouples were attached to the 
front and back surfaces of the mirror to measure the thermal contours produced. 
Figure 4-6 and 4-7 show the thermocouple distribution on the back and front 
surfaces of the mirror, respectively. The thermocouples were covered by small 
patches of pressure sensitive aluminum tape to help make good thermal contact 
with the mirror. Eight more thermocouples, which do not show in figure 4-6 
and 4-7 were spaced at equal intervals around the edge of the mirror. 
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Figure 4-5. Recommended Correlation for Free Convection 
Around Vertical Plane Surfaces. 

(From W.H. McAdams, Heat Transmission, 

3rd ed.. New York, McGraw-Hill, 1954.) 
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Figure 4-7. Thermocouple 
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on Front of 30 Inch Mirror 


On the back of the mirror, over the thermocouples, the heating 
blanket was cemented in place using contact cement. The arrangement and orien- 
tation of the heating zones and the thermocouple numbering arrangement are 
shown in figure 4-8. The guard strip which is divided into four equal heating 
zones was cemented in place over the thermocouples around the mirror edge. 
Thermistors were attached to each heating zone to monitor the individual zones 
for use in automatically maintaining stead temperature levels. 

Several layers of super insulation were placed between the heating 
blankets and the backing plate in order to reduce heat loss from the blanket 
to the backing plate and to isolate the mirror from the heating effects of the 
actuator motors. 

Twenty-six heaters were spatially arrayed on the mirror and indi- 
vidually held via associated thermiator sensors to desired set temperatures. 
Accurate determination of mirror temperature was accomplished by means of the 
array of thermocouples attached to the mirror front and back surfaces and its 
edge. Figure 4-9, a view of the associated control panel, shows the thermo- 
couple jack panel on the left. The right-hand half of the panel pictorially 
presents the shape and location of the individual heaters as they would appear 
when the rear of the mirror is viewed. Heater No. 1, for example is pie-shaped 
and located just below the center of the mirror. As is apparent from the figure, 
the heater locations are panel designated by numbers appearing on the white 
squares. The positions of the various thermocouples are designated by the 
other numbers. All numbers are on overaly sheets that can be lifted to permit 
access to the temperature set potentiometers, which are mounted to the panel. 

The dark circle below each heater designation number is the posi- 
tion of a neon indicator light, which lights whenever the blanket is receiving 
power. Directly below these neons are the screw-driver adjusted temperature 
set potentiometers. 
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Thermocouple Numbers 1 through 80 Are on Mirror Back Surface. 
Thermocouple Numbers 81 through 88 Are on Mirror Edge. 
Thermocouple Numbers 89 through 105 Are on Mirror Front Surface. 


Mirror Reaction 



Figure 4-r 8. Thermocouple Arrangement as Seen from Back of Mirror 
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Each of the 26 heaters was held at its set temperature by the 
control circuit as shown in figure 4-10. The arrangement features an A742 
I.C., controlling an L2001M7 four quadrant triac and operates in zero crossing 
switching mode to avoid unnecessary EMI generation. A resistor bridge, com- 
prised of the three 33K resistors plus an N.T.C. thermistor and series potentiom- 
eter, furnishes the necessary temperature error signal in electrical form. 

The separation between the two temperatures at which the control 
circuit turns on or off at a given setting, which is a peak to peak measure of 
how close the controller can maintain a given temperature at the thermistor, 
is 0 . 0 6 °F . 


In practice, the effectiveness of the control circuit was modified 
by the fact that the thermistors were placed on the back of the heating blanket 
instead of between the heating blanket and the mirror, where the thermocouples 
are. This was done in order to make better contact between the mirror and the 
blanket; however, it does mean that the temperature being regulated is the tem- 
perature on the back of the blanket. The effects of this were minimized to 
reduce heat loss, but the effects were still noticeable in trying to achieve 
the distributed profiles desired and were the cause of most of the non-uniform- 
ities observed in these profiles. 

4.2 TEST METHOD 

Information concerning the effects of local gradients was obtained 

as follows: 

(a) The mirror was aligned, the actuator positions and thermo- 
couple readings were recorded and the servo control was then turned off. 

(b) The heating elements were set for the desired gradient and 
the mirror was allowed to reach thermal equilibrium at the new gradient. 

(c) The mirror deformation was photographed and thermocouple 
readings recorded. 
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Fenwal Thermistor GA45P2 


Figure 4-10* Automatic Heater Temperature Control 


(d) The heating elements were turned off and the mirror allowed 
to reach stable equilibrium again at the initial thermal conditions. 

(e) Step (a) was repeated without turning off the servo control. 

(f) Step (b) was repeated with the servo control on to maintain 
alignment, and the actuator positions were recorded when thermal equilibrium 
had been reached. 

(g) The applied thermal gradients were obtained by subtracting 
the thermocouple readings in step (a) from the thermocouple readings in step 

(c). 


Information concerning the effects of distributed gradients was 
obtained in the following manner: 

(a) The mirror was aligned and the actuator and thermocouple 
position were recorded. 

(b) The heating elements were set for the desired gradient and 
the mirror was allowed to reach the desired thermal profile with the servo 
control on, to maintain alignment. 

(c) The actuator and thermocouple temperatures were recorded. 

(d) The servo control loops were opened and the heater elements 

turned off. 

(e) After overnight normalization, the resulting interferogram 
was photographed. 

4.3 RESULTS 

In practice it was difficult to approximate some of the distri- 
buted profiles shown in figure 4-1. In particular the projected difference 
between front and back profiles (the axial gradient) actually obtained was 
much smaller than expected in comparison with the radial gradients. 
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The actual temperature profiles obtained for the axisymmetric 
thermal distributions applied are shown in figure 4-ll(a). The measured data 
was averaged to give the symmetrical radial component for each of the four 
distributions. The temperature profile obtained for the non- symmetrical ther- 
mal distribution is shown in figure 4-12(a) . 

The temperature profiles obtained for the localized thermal gra- 
dients applied are shown in figure 4-13. 

The mirror deformations produced by the applied thermal gradients 
are obtained from the inter ferograms shown in figure 4-14 for the distributed 
gradients, and in figure 4-15 for the localized gradients. Figure 4-ll(b) 
shows the averaged symmetrical radial component for the four axially symmetri- 
cal distributed gradients. Figure 4-12 (b) shows the resultant deformation for 
the non- symmetrical distributed gradient of figure 4-12(a). Figure 4-16 shows 
a typical deformation profile for one of the localized external gradients . 

The changes in temperature measured at each of the thermocouple 
locations, from which the changes in thermal gradient were plotted, are given 
in Table 4-1 for the profiles shown in figures 4-ll(a) and 4-12 (a), and in 
Table 4-2 for the profiles shown in figure 4-13. The mirror figure changes 
induced by the changes in thermal gradient are measured from the inter ferograms 
and are given in Table 4-3 for the errors caused by the profiles of figures 
4-11 (a) and 4-12 (a) and in Table 4-4 for the errors caused by the profiles of 
figure 4-13. The forces applied by the actuators to remove the figure errors 
produced by the thermal profiles shown in figures 4-ll(a) and 4-12(a) are given 
in Table 4-5. 

The flex actuators had sufficient range to maintain the mirror 
figure for the gradients applied; however, the rate of heat application had to 
be controlled since the actuators could not always keep pace with the change 
when the heat was applied at the full 1 watt per square inch capability of 
the heating blanket. In some cases the maximum rate of heat application had 
to be reduced by a factor of two. 
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(a) Non-Ax lsyone trie Distributed Thermal Gradient 



Figure 4-12. Distributed Thermal Gradient Number 5 
and Resultant Mirror Deformation 
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(a) From Profile No. 1 


(b) From Profile No. 2 




(c) From Profile No. 3 


(d) From Profile No. 4 



(e) From Profile No. 5 


Figure 4-14. Interferograms Showing Mirror Deformations Due to 
Distributed Thermal Gradients of Figure 4-l(a) 
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(c) From Profile No. 3 (d) From Profile No. 4 


Figure 4-15. Interf erograms Showing Mirror Deformations Caused 
By Local Thermal Gradients of Figure 4-13 
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(a) Local Thermal Gradient No, 1 


to °C 

JL 



(b) Mirror Deformation Profiles 



Figure 4-16. Localized Thermal Gradient No. 1 and Resultant 
Mirror Deformation 
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TABLE 4-1 CHANGES IN TEMPERATURE MEASURED AT THERMOCOUPLES FOR THERMAL 
GRADIENT PROFILES SHOWN IN FIGURES 4-1 I (a) AND 4-12(a) 


Thermo- 

Couple 

Number 

AT ( ° C ) for Figure 
4-ll(a) Profiles 

AT for 
Figure 
4-12(a) 
profile 
#5 

Thermo- 

Couple 

Number 

AT(°C) 
4-11 (a) 

for Figure 
Profiles 

AT for 
Figure 
4- 12 ( a) 
Profile 

4 5 

4i 
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#3 
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TABLE 4-1. (Continued) 
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TABLE 4-2. CHANGES IN TEMPERATURE MEASURED AT THERMOCOUPLES FOR 
THERMAL GRADIENT PROFILES SHOWN IN FIGURE 4-13. 
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TABLE 4-2. (Continued) 
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TABLE 4-3. FIGURE ERRORS INDUCED BY APPLIED THERMAL PROFILES OF 

FIGURES 4-11 (a) AND 4-12(a) INFRINGES AT 6328& (1 FRINGE = 
12.46 MICRO INCHES). 

Displacement away from center of curvature defined as positive. 
(Actuator numbering arrangement as shown in Figure 3-5.) 
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TABLE 4-3. (Continued) 
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TABLE 4-4. FIGURE ERRORS INDUCED BY APPLIED THERMAL PROFILES OF FIGURE 4-13 IN 
FRINGES AT 63 28X (1 FRINGE = 12.46 MICROINCHES). 

Displacement away from center of curvature defined as positive. (Actuator 
numbering arrangement as shown in figure 3-5.) 
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TABLE 4-5. FORCES APPLIED BY ACTUATORS TO REMOVE FIGURE ERRORS PRODUCED BY THE 
THERMAL PROFILES SHOWN IN FIGURES 4-ll(a) AND 4-12(a). 

— O 

Forces given in pounds x 10 . (Actuator numbering arrangement as shown in 

Figure 3-5.) 
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TABLE 4-5. (Continued) 
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4.4 CONCLUSIONS 


The results show that the gradients produced by heat applied to 
small areas of the back of the thin deformable mirror are quite localized for 
pressures of one torr or less and the deformations produced by these localized 
gradients are also fairly localized. The effect of a distributed gradient 
appears to be equivalent to that which could be obtained by linear superposi- 
tion of the effects of a set of local gradients whose linear superposition adds 
up to the distributed gradient. Since, for the 30- inch diameter thin mirror, 
the effects of a local gradient produced by heat applied to the back of the 
mirror consistently displaces the corresponding area of the mirror surface 
away from the center of curvature, it should be possible to predict the effect 
of distributed gradients of low spatial frequency from the observed results. 

The deformation due to the distributed gradients appear consistent with the 
above conclusion. 
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5.0 


MODAL CONTROL 


5.1 INTRODUCTION 

It has been established that feedforward compensation can be 
used to decouple each of the control actuator channels while maintaining 
stability*. There are many possible approaches to optimization of the control 
system; however, a specific control law which takes advantage of the mirror's 
natural bending modes was selected for further investigation. 

The method of control, referred to as modal control, was sug- 
gested by J.F. Creedon and A.G. Lindgren** and they demonstrated the possibility 
of determining the minimum number of actuators placed in an optimum configura- 
tion so as to obtain a minimum rms displacement error, for a flat, freely 
supported plate subjected to a "white" spatial force distribution. Their re- 
sults show that 36 actuators could be employed for a 1/50 of a wavelength 

o 

(X = 0.6328A) rms error instead of 58 actuators as used in the initial phase 
of the thin deformable mirror project. 

One of the objects of the present phase of study was to determine 
similar information for the thin spherical mirror on a three point support. 

The bending mode shapes and their associated natural frequencies of vibration 
were analytically and experimentally determined (see appendix A) since this 
information is required to obtain the optimum actuator placement. These results 
yielded mode shapes and their associated frequencies for the three point support 
system. 


To determine actuator locations the strategy used is to control 
the dominant low frequency modes by placing control actuators at or near nodal 
contours of higher order modes which, if unexcited, do not contribute to the 
total rms error after control. Creedon and Lindgren found the square flat 
plate amenable to application of this strategy because the nodal contours fol- 

*Reference No* 1 
Reference No. 2 
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a relatively simple pattern. It was discovered from the resulting modal wave- 
shapes for the thin spherical mirror, however, that the nodal contours of the 
various modes did not follow a simple pattern when more than eight modes were 
considered. Further, it was also discovered that by including the 15 dominant 
modes, which comprised 99% of the specified mirror figure error of one wavelength 
rms amplitude, it is not possible at this time to determine their placement to 
obtain anywhere near \/50 error with control because of the very large number 
of possible actuator configurations. It became obvious at this point that a 
somewhat less ambitious approach would have to suffice to demonstrated that 
the modal control approach performed as predicted when applied to the actual 
active optics mirror and mounting system. 

The next phase of the study was to assume a limited number of 
actuator and sensor points to enable a design procedure using modal control. 

The approach taken was to select four control points and three additional un- 
controlled points. Thus, seven sensed points, including the actuator points, 
would feed data into ahmodal controller yielding four outputs which would drive 
their corresponding actuators with the desired result that: 

(a) The four controlled mode amplitudes should go 
from non-zero before control to zero after control 

(b) The three uncontrolled mode amplitudes should not 
change significantly. 

The former result would be obtained due to the presence of an 
integrating servomechanism if the actuators are not located exactly on the 
nodal contours of the first four modes. The latter result would only be ob- 
tained if the actuator placements were chosen at or near the nodal contours 
of the next higher order modes 5, 6, and 7. Experimental measurements verified 
that this approach gives the predicted results, thereby validating the modal 
control approach and furthermore establishing the design procedure as a workable 
one . 
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The next section deals with the details of the modal control 
approach in achieving the cited results. 

5.2 SYSTEM ANALYSIS 

A model of the system to be analyzed is shown in figure 5-1 of 

this report. A partitioned model is shown in figure 5-2. The model of figure 

5-2 will be used to show that the error vectors, e VT and e_, are identical to 
__ NR 

the mode vectors and C^. Then the simpler model of figure 5-1 will be used 
to calculate the mode vector, c, and the mirror figure vector, w, when modal 
control is being applied. 


From figure 5-2, the mirror acts as a mode analyzer which yields 
the mirror figure vector 



Likewise , 
tor , 


the controller contains a synthesizer which yields the control vec- 
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Figure 5-2. Partitioned Model o£ Modal Control System 
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Summarizing the analysis, the following resulting matrix equa- 


tions are: 
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(1) Controlled Modal Amplitude Vector: 

s ■ [i + % v s >r % • <w 

where q is the disturbance vector for the N 
N 

controlled modes. 

(2) Uncontrolled Modal Amplitude Vector: 

' % ' X ^ <S> ^ i' 1 [l + U ^C)] -1 q N i (EX) 

where q R is the disturbance vector for the R 
uncontrolled modes. 


(3) Actuator Displacement Amplitude Vector: 

"» = V» + 5* s 

(4) Uncontrolled Displacement Amplitude Vector: 


W„ 


U RN C N + U RR 


where U-_ T , U x _ , U_ XT , and U„ are partitions of the modal matrix U__; N repre- 
NN NR RN RR M 

sents the number of actuators being controlled out of M possible points and 
R * M - N. 


gain k. 


The above results assume simple integrating actuators of equal 


If the disturbance vectors q~ and q are step functions, the 

« R 

resulting steady state mode amplitude vectors can be calculated from (1) thru 
(4) as: 


(5) 



= 0 


( 6 ) 



’r • v°> v v " 1 v°> % 
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<7 > \ s - V C K„ S 


< 8 > W R * »u °R 

ss ss 


The rtns mode amplitude in steady state will then be 


(9) C 


RMS 


» ixf 


ss 


(M = the number of sensor points) 


Likewise, the rtns figure error will be 


( 10 ) 


“rhs - 5 X V J »ss ■ 5 / (< V + V )C R 

I * ” ss ' 1 1 SS 


5.3 VERIFICATION OF MODAL CONTROL APPROACH 

Having established the validity of the mirror's natural modes 
and characteristic frequencies it was decided to perform a design for a test 
case on a limited number of actuators, with experimental verification. The 
work performed herein gives the analysis leading to the design of a four- 
actuator system. Seven sensing points that include the four actuator points 
were chosen. It is the design objective to diminish the four most dominant 
modes by modal control, without simultaneously exciting the three next most 
dominant modes. 

5.3.1 Selected Configuration 

Figure 5-3 shows sketches of the nodal contours (locus of zero 
deformation) of the seven most dominant modes, 1 through 7. As the mode num- 
ber increases, characteristic frequency increases and the corresponding ampli- 
tude tends to decrease monotonically . Figure 5-4 shows a composite or super- 
position of modes 5, 6, and 7, which are not to be excited while controlling 
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Mode 1 



Figure 5-3. Nodal Contours of First Seven Dominant Modes 
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Mode 5 



Figure 5-4. Superposition of Modal Contours for Modes 5, 6, and 7 



modes 1 through 4. It is then obvious from figure 5-4 that the circled points 
are the actuator locations. The three cross-marked points are chosen as 
additional sensing points. Thus, the actuator points chosen are 20, 17, 23, 
and 47 according to the numbering scheme shown in figure 5-5 and three addi- 
tional sensor points are chosen at 6, 41, and 37. 


5.3.2 Design Procedure 

With the actuator and sensor points chosen it was then possible 

to determine the modal matrix U„, (7x7 since N = 7) which is a submatrix of 

N 

U (58x58 since M = 58) given as 
M 


u 


M 


u^l) 

u 2 (l) .. 

'•• u 58 (1) 

u 1 (2) 

u 2 (2) .. 

...u 58 (2) 


Ul (58) u 2 (b8) u 58 ( 58 > 


( 1 ) 


where the notation u.(i) refers to conventional u. . notation. That is, the 

3 ij 

number is in parentheses refers to the i-th point in question and the subscript 
number j refers to the j-th mode. Row 1 in equation (1) is column 1 in the 
conventional sense, while column 1 in equation (1) is row 1 in the conventional 
aensev. Thus, for example, u^ " u 2 q 2 = u 2 ^^) since the chosen point 1 for 
the experiment corresponds to grid point 20, which is the center point of the 
mirror. 


Written out in detail the modal matrix, U w , is 
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The dotted lines refer to partitions that define, respectively, 
from top to bottom or left to right, actuator points* Comparing equation (1) 
and equation (2) yields a relation between the natural order in equation (1) 
and the actual order in equation (2) as shown in Table 5-1. 


TABLE 5-1. POINT CORRESPONDENCE 


Natural Order 


Actual Order 


1 

2 

3 

4 

5 

6 

7 


20 

17 

23 

47 

6 

41 

37 


etc. 


etc. 
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The displacement vector w ia related to the mode amplitude vector c by 


(a) 


V m U c or 


( 3 ) 


(b) c»U _ 1 w«Vw ; v - u 1 (58 x 58) 


In accordance with the partition mode 


"N 


V NN W N + V NS W s 


C - V CM W M + V C o w o 

s SN N SS s 


(4) 


where C„ and C refer to the mode amplitude vectors for the actuator points 
Ns _ _ 

and sensor (non-actuator) points, respectively* and W refer to the actuator 

H s 

and sensor points, respectively. The V partition matrices are obtained from 

V = U -1 . By referring to Table 5*1 this establishes the submatrix V _ as 

N+S 

(5) 
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From numerical data. 



I 

-4.3155E-3 

-7.9662E-1 

-8.0690E-2 

8. 6487E-1 

(a) 

V * 

NN 

2.1392E-3 
-5 . 7397E-1 

-5.4387E-1 

-5.4445E-1 

9.6133E-1 

-5.4421E-1 

-4.1127E-1 

-5.5683E-1 


! 

-7.8007E-1 

4 . 6961E-1 

4. 6960E-1 

4 . 6991E-1 



"-1.8454E-1 

1 . 6137E-1 

1 . 2446E-2 " 


(b) 

V NS = 

8.7841E-2 

-4.5913E-1 

1.1592E-1 

-4.6118E-1 

-1 . 9782E-1 
-4. 6122E-1 




. -4.2135E-1 

-4.2454E-1 

-4.2460E-1 j 




" -1 .4855E-3 

-4.3910E-4 

-4.3316E-4 

6.8879E-3 

(c) 

V SN ' 

-1 . 1527E-5 

-1.2895E-3 

-1.2521E-3 

-9.0923E-6 



-2.5437E-5 

4.8284E-3 

-4.7930E-3 

2.4387E-5 



-7 . 7551E-1 

3.8640E-1 

3.8551E-1 


(d) 

v ss 

1 . 61898E-3 

-1.8611E-1 

1.84474E-1 




-1.02063E-3 

-6.4466E-1 

6.4564E-1 



The physical realization of a coefficient set can be obtained 
by the op-amp resistor network configuration shown in figure 5-6. From figure 
5-6, the output voltage obtained is 
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Figure 5-6. OP-AMP Resistor Network Realization of Controller 


where the minus and plus superscripts refer to the negative and positive coef- 
ficients, respectively. The displacement inputs w^ really are proportional 
phase detector output voltages. Since. 


then 


or 





M 

£ 


a u w j 


( 8 ) 


(9) 

( 10 ) 


thus, for a suitably scaled feedback resistor, Rp, the input resistor values 
are essentially proportional to the reciprocal of coefficients of the V matrices 
of equation (6). The resulting mode analyzer realization is shown in figures 
5-7 and 5-8. 
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Report No, 10859 


Finally, the effect of actuator pads* is considered by intro- 
T 

ducing the pad matrix, H « GU . For small pads G = 1, the identity matrix. 

T 

Assuming this to be the case, H * U . For pad compensation, the required con- 
troller must be \ where the subscript N refers to the actuator points only. 
This matrix was calculated and resulted in the network of figure 5-9. 

5.3.3 implementation of Controller Design 

The actual implementation of the controller differed only 
slightly from that shown in figures 5-7, 5-8, and 5-9. Computation of thru 
did not require the use of the six inverting amplifiers as indicated, since 
the phase detector output signals as well as their inverted counterparts were 
available from existing electronics. Associated with the C^, 5 and 

generating operational amplifiers k^, ky and A^) are four inverting 

amplifiers (A A A_^ , and that generate -C^, -C ^ > -C 3 , and -C^. 

Hence, the generation of subsequent signals (via pad compensation) was accom- 
plished without the need for "inverting" amplifiers as shown in figure 5 - 9 . 
However, a summing amplifier was required in each of the four output channels 
to drive the associated servo amplifier summing point through a 10K output 
resistor. This added complexity was required to enable use of the large value 
summing resistors shown in figure 5-2. This prevents overloading of the pre- 
ceding amplifiers as might otherwise occur. 

Figure 5-10, a half size phusical layout drawing, shows the 
implemented circuitry. It should be noted that resistor values up to 140K 
are implemented by the nearest 1% resistor value. Higher resistances between 
1 and 10 megohms are implemented with a 50K potentiometer whose slider feeds 
a portion, 1 /x, of the applied voltage signal, ±W^ , to a one megohm summing 
resistor. This implements an equivalent summing resistor value of X megohms. 
The same approach is used to obtain the equivalent of summing resistors between 
10 and 100 megohms; however, here a 10 megohm fixed resistor is used with a 
500K potentiometer. Even higher equivalent summing resistors (>100M ) are 

*J.F. Creedon and A.G. Lindgren, Op. Cit. 
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obtained in essentially the same manner. There a 50K potentiometer feeds a 
100:1 divider whose output voltage drives a one megohm fixed resistor. Cermet 
potentiometers capable of 1% adjustment are used for compatibility with the 1% 
fixed resistor covering the <140K range. 

It should be noted that the small circles labelled and ±C^ 
represent possible locations for the potentiometers associated with summing 
resistors. A value (either 50K^ or 500 K^) adjacent to the location indicates 
that a (50 or 500K ) potentiometer was included to drive the resistor located 
to the right. In cases where no value is shown adjacent to the small circle 
location, a potentiometer mounting hole but no potentiometer was provided. In 
these cases, the input signals were applied directly to the summing resistor 
to the right. All amplifiers except those used for signal inversion feature 
two zener diodes paralleling the feedback resistor, to provide clipping action 
in all cases at about 9 volts. 

The controller as finally breadboarded for experimental testing 
is shown in figure 5-11. The input and output wires were cabled to connectors 
that mate with the existing connectors at the phase detector outputs and servo 
amplifier inputs so that the unit can be inserted into or removed from the 
system quickly and conveniently. 
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Figure 5-11. Four Actuator Modal Controller Breadboard Unit 


5.3.4 Experimental Results with Modal Controller 


The modal controller was tested by introducing errors into the 
previously aligned, 30-inch mirror and then removing as much as possible of 
the error introduced, with the modal controller. The relative amplitudes of 
the modal error signals produced by the uncontrolled error were observed and 
compared with the modal amplitudes predicted analytically for the specific 
error Introduced. The modal amplitudes remaining after closed loop correction 
were measured. The time responses of the system were observed and the change 
in forces required to effect correction were recorded. The remanent error 
after correction was also compared with the remanent; error produced by the 
same four actuators in parallel closed loop operation without the modal con- 
troller. 


The results of the measurements are described in the following 
pages for two different error disturbances: one, a point force error introduced 
by a single actuator applying a force against the mirror to distort it; and 
two, a spherical error introduced by changing the axial distance between the 
mirror and the figure sensor. 

The point force was exerted at actuator number 28 (see figure 
5-5). This is an uncontrolled point halfway between a controlled point (actua- 
tor number 2Q) and sensor point (actuator number 37). 

The errors developed as a function of the force applied, without 
correction, are shown in figure 5-12. The phase detector outputs are numbered 
through corresponding to the seven sensor points. The modal controller 
outputs, which should be proportional to the amplitude of the first seven 
natural modes contained in the mirror deformation are numbered through C^. 
The relative slopes of through are compared in Table 5-2 for the measured 
values and the predicted values and show good correlation. 
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Open Loop Phase and Mode Error Signals (Volts) 



Figure 5-12 • Open Loop Phase and Mode Errors as a Function of 
Force Applied at Uncontrolled Actuator No* 28 
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TABLE 5-2 


C 1 

Measured 

1.0 

Calculated 

1.0 

C 2 

-10.5 

-16.4 

C 3 

-52.0 

-52.6 

c 

-10.0 

- 6.12 

4 

C 

2.0 

1.73 

5 

C 6 

1.5 

0.83 

c. 

5.2 

2.89 


The phase detector outputs and the modal controller outputs 
after closed loop correction for the point force at actuator 28 are shown in 
figure 5-13. Since the phase detector outputs are only linear up to errors of 
1/8 wavelength, the system can only be expected to operate properly within 
this limit. It is interesting to note, however, as shown in figure 5-13, that 
the system operates stably, if not accurately, for errors greater than 1/8 
wavelength up to 3/8 wavelength amplitude. As expected, modal controller out- 
puts through are nulled identically. The change in forces generated in 
the four actuators to maintain alignment are also shown in figure 5-13. 

The mirror deformation introduced without closed loop control 
for an applied force of 0.175 pounds at actuator number 28 is shown in figure 
5- 14(a). Figure 5- 14(b) shows the residual error after nulling through 
with the modal controller. For comparison, figure 5-14 (c) shows the residual 
error after nulling the four actuators using the straight parallel control 
loops. 


The errors developed without active correction as a function of 
change in axial spacing between figure sensor and mirror are shown in figure 
5-15. 
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Figure 5-14. Mirror Deformation Produced by 0.175 lb. Force 
Applied by Actuator No. 28. 
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Open Loop Phase and Mode Error Signals (Volts) 
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Mode Error Signals 
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Change in Axial Spacing (Inches) 

Figure 5-15. Open Loop Phase and Mode Errors as a Function 

of Axial Spacing Between Mirror and Figure Sensor* 
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The phase detector outputs, the modal controller outputs, and 
the actuator forces developed after closed loop correction are shown in figure 
5-16. 


The mirror deformation introduced without closed loop correction 
for a change of 0.003-inch in axial spacing between the mirror and figure sen- 
sor is shown in figure 5-17 (a). Since the error is axially symmetrical, a 
single diameter profile is shown. Figure 5-17 (b) shows the residual error with 
modal control and for comparison figure 5-17 (c) shows the residual error for 
the straight parallel closed loop control. 

Both cases shown exhibit improvement in RMS figure for the 
modal controller with four control actuators and seven sensing points, over the 
straight parallel control with four actuators and four sensing points. The 
larger the number of sensing points used, the better the fit should be to the 
reference wavefront, for any given number of actuators. 
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3 


(a) Without Closed Loop Control 



(b) Modal Control with Four Aotuator* 



(c) Straight <-r*xa 11*1 Oo*tr*l wick Four dot oat or a 

Figure 5-17* Mirror Figure Error Produced by 0.003” Change in 
Axial Spacing Between Mirror and Figure Sensor 


6.0 TRANSIENT RESPONSE TO STEP ERRORS 

Tests were made on the three different control configurations; 
the straight-parallel coupled control loops, the uncoupled feedforward network, 
and the modal controller, to determine their transient response to step errors. 
The step errors were obtained by changing the axial spacing between the mirror 
and figure sensor or by displacing the individual flex actuators by various 
amounts. When possible, the control loops were first opened, the step dis- 
placement errors were created and the control loops were then closed. The 
responses were measured and recorded at the output of the phase detectors from 
the time of closure of the loops, except for the mode errors which were read 
at the modal controller operational amplifier outputs. 

Response curves are shown in Figure 6-1 for the straight-parallel 
coupled control loops. Figure 6-2 shows the response with the feedforward net- 
work, and Figure 6-3 shows the response with the modal controller. 

The straight-parallel control loops exhibit a uniform response 
to an axial step error in contrast to the differences observed between the 
responses of actuators at different locations when the error was due to indi- 
vidual actuator displacement errors. The time constant for the axipl step is 
much shorter than the time constant for response to individual actuator dis- 
placement errors because all actuators see an error of the same sign and move 
together in the same direction to remove it, so that they are helping each 
other, giving an apparent increase in loop gain. 

Parallel control loop configurations using fewer actuators 
spaced further apart also give shorter time constants because of higher effec- 
tive loop gains. 

In contrast, the response of the feedforward network is not 
only uniform from one actuator to another for individual actuator errors, but 
it is essentially the same for axial errors. This is because the gain has been 
normalized and decoupled and the actuators are in effect unaware of whether 
the adjacent actuators are aiding or opposing them. 
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58 Actuators - Local Step at No. 20 



58 Actuators - Local Step at No. 16 



58 Actuators - Axial Step (Typical) 



4 Actuators - Axial Step (Typical) 



Figure 4 - 1 . Actuator Response to Step Errors with Straight- 
Parallel Coupled Loops 
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Figure 6*2 # Typical Actuator Response to Axial Step 

Error or Individual Actuator Displacement 
with Feedforwards. 
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7.0 LOCALIZATION OF FEEDFORWARDS 


The feedforward network is a resistive matrix designed to com- 
pensate for interactions in the deformable mirror so that the individual con- 
trol loops can operate independently. The design of this network is described 
in detail in Reference 1. The objective here was to determine how well the 
design values come to giving the desired results. The feedforward network 
is designed to give a localized deflection on the mirror if an error signal is 
applied to a single input of the feedforward matrix, so that a measure of the 
network's effectiveness is its ability to produce a localized displacement. 

The test for the ability to produce a local displacement was to put a DC volt- 
age at the input of the feedforward matrix, one channel at a time, with the 
control loops open, and to observe the resultant deflection as indicated by 
the interferometer output. Some typical interferograms are shown in figure 
7-1. Figures 7-l(a) and (b) show the desired result for a signal applied to 
control loop number 40. The deflection actually produced with the network as 
originally designed, however, is shown in 7-l(c) and (d) . The effect shown 
in 7-1 (a) and (b) was obtained by decreasing the value of driving point resis- 
tance approximately 30 percent from the design value. The reason this was nec- 
essary was discovered by observing the force changes in the actuators involved. 
The second ring of actuators stayed essentially in their original positions even 
though a six volt signal was applied to the driving point. This means that the 
full set of 19 forces designed to produce local displacement was not being 
applied. This was caused by the actuator dead zones which, for the actuators 
in the second ring, are equivalent to approximately ±1/200 wavelength. The 
driving point resistance correction made to obtain 7-l(a) and (b) gives essen- 
tially a single ring (7 actuators) constellation of forces for each control 
loop instead of the 19 point group originally designed. The localization ob- 
tained with this more limited compensation is, however, very good as can be 
seen in figure 7-l(a) and (b) , and in the light of hindsight it is evident that 
the 10 point arrangement was over compensation even if the dead zone phenomenon 
did not exist. 

The dead zones are approximately twice as large with the feed- 
forward network as with the straight parallel control because the minimum driv- 
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(C) 


w 


Figure 7-1 » Interf erograms Showing Localized Deflections Produced By 
Feedforward Network (a) and (b) After Decreasing Driving 
Point Resistance to Produce a Seven Point Configuration 
of Forces (c) and (d) With Original Design Values for 19 
Point Force Configuration 
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ing point resistances are at least twice as large. The dead zones for the 
straight parallel control loops are equivalent to ±1/400 wavelength or less. 

Some dead zone is usually desireable to prevent overworking the actuator motors. 
If friction did not provide a dead zone it might be necessary to provide a 
threshold control. 

An attempt was made to measure the stiffness matrix values by 
displacing a single actuator and measuring the change in actuator forces re- 
quired to maintain all the other points aligned. A typical set of forces is 
shown in figure 3-2. The dead zones here again make it difficult to obtain a 
set of readings that would be accurate enough to improve on the analytical 
determination. 
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APPENDIX A - DERIVATION OF MIRROR SUPPORT SYSTEM El GENFUNCT IONS 


In order to develop a controller design using the modal control approach^ 
it was necessary first to determine the eigenfunctions for the spherical mirror 
and its three point support system. This was accomplished by three different 
methods which are: 

1. Analytical solution of the beam equation for a thin flat circular 
plate subject to boundary conditions. This is an approximate solution since 
the actual mirror is a thin shallow shell. 

2. Finite difference approximation to the continuous medium and solution 
obtained using the SAMIS computer program. 

3. Experimental verification of 1 and/or 2 using time-lapse holography 
upon the physical system set into vibration. 

The results of these three studies was presented recently at the April, 1971 
meeting of the Optical Society in Tucson, Arizona. The presentation given, which 
is contained herein, establishes excellent agreement between the three approaches. 
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A. 1 ANALYTICAL APPROACH 


There are four types of vibrational inodes when a homogeneous, 

isotropic, thin circular plate is simply supported at three points on the 
* 

circumference. 

The type I is, 
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where k. is the i^ solution from 
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The type II is. 
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where k^ is the i-th solution from 


(A-2) 


00 

1 a (ka) 
m-l, 2,4, 5,7,8... m 
m^O.3,6,9,12 . . . 


- 0 


The type III is. 
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*The derivation of these results can be found in Perkin-Elmer Report No. 
10299, Evaluation of the Thin, Deformable Active Optics Mirror Concept , 
No. 2, Sept. 15, 1970. 
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where is the i-th solution from 
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The type IV is. 


(A-4) 


w - [J (k p) + D I (k p) ] sin mcp 

Ivx # m ml m m l 

where m - 3,6,9,12... (A-5) 

and where k is i-th solution from 

F (ka) - 0 

m 

The resulting mode shapes are shown in Figures A-l to A- 12 


One notes that types I and II are the even modes and types III 
and IV are the odd modes with respect to the x-axis. One notes also that the 
types I, II and III exert forces on the support points, whereas the type IV 
does not exert forces on the support points. Therefore, the types I, II and 
III require the support points in order for them to exist, and they are the 
modes which are unique to the plate with the simple support points on the 
circumference. The type IV is the group of modes which are also members 
of vibrational modes of the free plate without support points. These modes 
do not require the support points in order for them to exist. 
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Figure A-l Deflection Mode, Type I, k = 2.385, Even Mode 
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Figure A-2 Deflection Mode, Type I, k =* 3.289, Even Mode 
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Figure A-4 Deflection Mode, Type II, k » 2.453, Even Mode 
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Figure a- 8 Deflection Mode, Type III, k = 2.453, Odd Mode 
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A. 2 EXPERIMENTAL DETERMINATION OF MODAL CONTOURS BY TIME-LAPSE HOLOGRAPHY 

The method of time-lapse holography was used to experimentally 
determine information about the thirty-inch mirror's vibration mode shapes. 
Ideally, each mode shape occurs at a characteristic frequency. 

A time-lapse holograph is a film record of the combination of a 
reference beam and a time varying object beam. The principle involved is 
that the time average of a vibrating object's wavefronts, which are emitted 
continuously, is exactly equal to the average of the sequence of the object's 
wavefronts taken simultaneously. The validity of the latter statement is clear 
if we assume that an ordered sequence can be found from the continuous emission 
of wavefronts. 

In a conventional test setup, a hologram is made of the object in 
a reference state. The developed hologram is then inserted into its original 
position when first taken. The object is then deformed so as to produce 
optical interference between the object's deformed wavefronts and its station- 
ary wavefronts. Upon illumination by both the reference and object beams, an 
observer will see the reconstructed image with interference fringes representing 
the deformation superimposed on the image. If the reference beam is blocked, 
only the interference fringes that are due to deformation will appear. If 
the object beam is blocked, the reference beam will cause only the stationary 
object to appear. 

This approach, somewhat modified, was taken to obtain interference 
patterns concerning the mirror's vibration modes. The modification is that 
the step of making a reference hologram was deleted since it was not necessary 
when using active optics. The automatically aligned object mirror, together 
with the Twyman-Green interferometer portion of the figure sensor, establishes 
the reference wavefronts. The interference of these reference wavefronts 
produces a uniform field that holds no information. Thus, upon exciting 
a mirror mode at its characteristic frequency, the resulting interference 
pattern contains the information necessary to determine the mode shape. A 


A- 16 



contour plot of the mode shape is taken from the interference pattern because 
of the relationship between displacement amplitude and recorded light intensity. 
This relationship is fundamental in all fringe analyses, and its derivation 
is included in the appendix. The relationship is that the film exposure 
(intensity x time) E, is related to the amplitude of a sinusoidal vibration 
A, by 


E 


E 

o 


( 


1 + mJ 




(A-6) 


where E q (1 + m) is the average exposure with no vibration (A-0) ; m is a non- 
negative modulation factor less than unity (see appendix) ; and is the zero 
order Bessel function. A sketch of the exposure as a function of vibration 
amplitude is shown in Figure A- 13. 


E 



Figure A- 13 Normalized Exposure as a Function of Vibration 
Amplitude for a Time Lapse Hologram 
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Figure A-14 shows a positive interferogram taken at an excitation 
frequency of 272 Hz. This picture verifies the exposure relationship given 
by Eq.(A-6)in Figure A-13 

It is first noted that the brightest triangular contour corresponds 
to the three known support points of zero axial displacement . As we move 
radially inward, the fringes appear with diminishing brightness, thereby 
indicating monotonically increasing amplitude. This fact is explained 
physically by considering that points of motion move from the stationary, 
bright, uniform field to a dark field and then back to light, and so on. 

Upon averaging light intensity at each point, a gray level will result 
depending upon the time spent in bright fringes and dark fringes. Points 
of little motion spend most of their time in the original bright field, 
while points of increasing displacement average out to decreasing modulation. 

From analysis of Figure A-14 } a contour plot of the mode shape 
was made and is shown In Figure A-15. To test the validity of this mode shape, 
a comparison was made with preliminary results obtained from the SAMIS 
computer program. The comparison shows good agreement with SAMIS mode 04 
(radial-azimuthal) sketches in Figure A- 16 

Other Interesting interf erograms were made at other characteristic 
frequencies; some are shown in Figures A-17 f A-18, and A-19 (More comparisons 
and checks will be made for these other modes to ascertain the validity of this 
data. ) 


The word Interferogram is used here interchangeably with hologram since an 
interferogram is a special case (on-axie parallel beam) of a hologram. 
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Figure A- 14 Interferogram of Vibrating Mirror at 272 Hz 
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Figure A-17. Interf erogram of Vibrating Mirror at 296 Hz 
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Figure A- 18 Interferograms of Vibrating Mirror 



A»3 RELATION BETWEEN MODAL CONTROL LAW AND MIRROR STIFFNESS CONTROL LAW 


This section will show that when the number and placement of con- 
trolled displacement points and sensing points are identical, the controller 
design for perfect decoupling is the same whether the modal control approach is 
used or whether the original control law is used. The original control law 
used the mirror system* s stiffness matrix to cancel the displacement interaction. 
The difference between the two systems is that the modal control law takes 
into account displacement error between the controlled points while driving 
the lower mode amplitude to zero; the stiffness control system drives displace- 
ments at the controller to zero without regard to displacement error between 
the control points. 

To show the relation between the two methods let us consider the modal 
control approach. Following Creedon's notation we may expand the figure's 
displacement, w, and applied loading, p , according to 


w (x,y t) = c i (t) x U i (x,y) 

P (x,y,t) = St a. (t)^ ( x ,y) 


( 1 ) 


where ^/Wj(x, y) ^ is the orthonormal set of eigenfunctions of the mirror 
and support system. Their corresponding eigenvalues, “ 7 \ ^ , relate the mode 
amplitudes by 

c. = 'x . a i 

i ''x i 

(i = 1, 2, 3, . ♦ .m) 


( 2 ) 
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or in matrix notation 


C = A a (3) 

where c and a are the mode amplitude column vectors of dimension M x 1 = 58 x 1 
and A I s the eigenvalue diagonal matrix of dimension M x M = 58 x 58, Writing 
Eq. (1) in matrix notation 

W = U c (M x 1 = 58 x 1) 

(4) 

P = U a (M x 1 = 58 x 1) 

where each column of U (dimension M x M ° 58 x 58) represents an eigenfunction 
whose components are a mode amplitude evaluated at the sensor points. 

From equations (3) and (4) : 

C = _AI 1 p (5) 

It is also known that the influence coefficient matrix, a (dimension M x M = 

58 x 58), relates applied load to resulting displacement by 
W = a P = U C (by Eq. (4)) 

Then 

c = u" 1 a p (6) 

Equating Eq. (5) to Eq. (6) 

u -1 a = A u" 1 (7) 

or 

a = u A u" 1 
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But since the eigenfunctions form an orthonormal set (not complete) 

if 1 = U t (transpose) (8) 

so that 

o( = U V 1 (9) 

Equation (9) is the key relation between the two control methods. 

In the stiffness decoupling approach the controller is represented by the stiff- 
ness matrix 


k = or 1 = (if )” 1 A ~ l u -1 = u A -y 1 

But the modal controller according to figure 5-1 is 


( 10 ) 


- “ n " 1 “(S )-” 1 


( 11 ) 


where is the pad compensation matrix which prepares the load distribution 
to be a scalar multiple of a resulting displacement for any specific mode. 
Since H N = U t for small pads (see reference 2), 

V 1 = (uV 1 = U, 

So that (12) 

£ - £ £(s) - 2 d < b) Jl‘ 


Thus , 
matrix 


the modal controller matrix, £ 


9 ]£ if — ** Other wise 


, is identical to the stiffness controller 
they differ by a scaler matrix. 
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